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ABSTRACT 

We show, theoretically and via MHD simulations, how a short burst 
of reconnection localized in three dimensions on a one-dimensional cur- 
rent sheet creates a pair of reconnected flux tubes. We focus on the 
post-reconnection evolution of these flux tubes, studying their velocities 
and shapes. We find that slow-mode shocks propagate along these recon- 
nected flux tubes, releasing magnetic energy as in steady-state Petschek 
reconnection. The geometry of these three-dimensional shocks, however, 
differs dramatically from the classical two-dimensional geometry. They 
propagate along the flux tube legs in four isolated fronts, whereas in the 
two-dimensional Petschek model, they form a continuous, stationary pair 
of V-shaped fronts. 

We find that the cross sections of these reconnected flux tubes appear 
as teardrop shaped bundles of flux propagating away from the reconnec- 
tion site. Based on this, we argue that the descending coronal voids seen 
by Yohkoh SXT, LASCO, and TRACE are reconnected flux tubes de- 
scending from a flare site in the high corona, for example after a coronal 
mass ejection. In this model, these flux tubes would then settle into equi- 
librium in the low corona, forming an arcade of post-flare coronal loops. 

Subject headings: Magnetic reconnection — MHD — Sun: flares — Sun: CMEs 



1. INTRODUCTION 

Coronal mass ejections (CMEs), sudden eruptions of coronal plasma and mag- 
netic field into the interplanetary medium, are one of the most dynamic phenomena in 
the solar corona. When these ejections collide with the Earth's magnetosphere, they 
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can have dramatic effects on the Earth's space weather. Even if they miss the Earth, 
these CMEs are an important driver of space weather through the magnetic recon- 
nection and flaring that is driven in the corona in their wake. Particles accelerated 
during the flares associated with CMEs can damage satellites and create radiation 
hazards for astronauts, and significant EUV and X-ray emission is generated by the 
coronal plasma heated by these flares. 

Theories for the cause and form of CMEs vary (see, e.g., reviews by Forbes 
2000; Klimchuk 2001; Low 2001; Lin et al. 2003), but in all cases the CME is a 
magnetically dominated structure, and magnetic energy release plays a prominent 
role in its eruption. In these models, magnetic field lines, carrying coronal plasma 
with them, bow out from the low corona to form the CME, as in the simulation shown 
in Figure 1 (MacNeice et al. 2004). 

A critical step of the process is the reconnection of magnetic field in the wake of 
the CME. Prior to such reconnection the inward and outward directed field fines in 
the two legs of the CME are pinched together to form a current sheet. When the field 
reconnects across this sheet, as shown in the right hand panel of Figure 1, the CME 
is untethered from the Sun while the field on the sunward end of the sheet forms 
post-eruption coronal loops. As more of the field reconnects, these post-eruption 
coronal loops build up on top of each other to higher and higher altitudes, and the 
footpoints of each new set of loops arc more and more widely separated. This classic 
two-dimensional flare model (sec, e.g., Carmichael 1964; Sturrock 1966; Hirayama 
1974; Kopp & Pneuman 1976) explains both the appearance of progressively higher 
coronal loops in a post-CME reconnection event, and the appearance of progressively 
more widely separated Ha footpoints of the loops. 

The basic process of reconnection across a current sheet was flrst studied in 
two-dimensional, steady-state models which assumed uniform resistivity (Sweet 1958; 
Parker 1957). The reconnection rates in such models were usually too low to explain 
the fairly rapid development observed in CMEs (see, e.g., Biskamp 1986). A faster 
reconnection regime was proposed first by Petschek (1964); it was distinguished by a 
very small diffTision region with four slow-mode shocks linked to it. It turns out such a 
structure is inconsistent with spatially uniform diffusion (Kulsrud 2001), but occurs 
naturally whenever the magnetic diffusion coefficient, or any other non-ideal fleld 
line transport mechanism, is locally enhanced. This explains the notable successes 
of a number of modified reconnection models whereby fast reconnection can occur 
across a current sheet provided it is spatially localized, as by a locally enhanced 
resistivity (Ugai & Tsuda 1977; Scholer & Roth 1987; Erkaev et al. 2000). The degree 
of enhancement is not nearly as important as its localization (Biskamp & Schwarz 
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2001). 

While the hypothesis of locahzation solves the puzzle of fast reconnection, it 
raises new questions about how such a small-scale phenomenon might couple to the 
global magnetic field, say of the CME. A partial answer was provided by investigations 
of two-dimensional, localized, non-steady reconnection beginning with Semenov et al. 
(1983) and subsequently developed by Biernat et al. (1987), Heyn & Semenov (1996) 
and Nitta et al. (2001). These models postulate a localized, non-steady reconnection 
electric field on a pre-existing current sheet and solve for the external response. The 
timc-dcpcndcnt solution, shown in Figure 2, consists of a pair of post-rcconncction 
loops retracting at the Alfvcn speed. Some models also include a small- amplitude 
fast mode wave propagating radially outward, in advance of the slow-mode shocks 
(Semenov et al. 1983; Heyn & Semenov 1996; Nitta et al. 2002); others such as Biernat 
et al. (1987) treat this as a static disturbance to the up-stream field. 

This two-dimensional model demonstrates several respects in which locahzed, 
non-steady reconnection differs from the more well-known steady-state scenario. For 
example, during reconnection four slow-mode shocks extend from the reconnection 
site as in the steady-state Petschek model, but they do not extend to infinity, as in 

the steady-state. Rather they close around the tips of the unrcconncctcd flux tubes 
(see Fig 2). In addition, when reconnection ceases the flux-tubes continue to retract, 
with the teardrop-shaped slow-mode shock fronts remaining intact, and a new current 
sheet forming behind them. Even then, the magnetic energy of the system continues 
to decrease as the flux tubes retract. This means that the energy released by the 
reconnection (i.e., converted from magnetic to kinetic and thermal energy by the 
shocks) can far exceed the energy actually dissipated (converted directly to heat by 
the reconnection electric field). Finally, the loops continue to sweep up mass into 
their slow-mode shock envelope as they retract, thereby growing indefinitely through 
snow-plow-like effect (Semenov et al. 1998). 

The present work is aimed at generalizing this scenario to localized, non-steady 
reconnection in three dimensions, which we call patchy reconnection. The objective is 
to characterize the effect on the global flcid of reconnection occurring within a small 
region of a pre-existing current sheet. In particular, we solve for the dynamics of 
the post-reconnection ffux tubes once they have been created by some reconnection 
process. This is a question of fundamental importance in the study of magnetic 
reconnection, as well as a critical element in understanding CMEs. 

Observational approaches have been slow to yield this understanding since post- 
reconnection flelds are not observable immediately after the CME. Rather, only the 
resulting changes in the coronal configuration are observed. This is because newly 
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forming reconnected loops only become visible in emission after most of their energy 
has been released and they fill up with hot plasma. This stage only occurs as the loops 
approach their equilibrium state in the low corona, and so both the initial reconnec- 
tion and the majority of the subsequent dynamics are unobserved. Thus even truly 
dynamic studies such as those of Forbes & Acton (1996) can only observe the end 
stage of the loop evolution. The apparent dynamics which the two-dimensional mod- 
els explain so well are simply the gradual brightening of successively higher coronal 
loops with successively wider footpoints. These two-dimensional reconnection models 
can therefore only compare their predictions with observations of the final equilib- 
rium of the magnetic field. While they have been quite successful in this endeavor, 
the dynamics of the reconnection and even of the loops themselves remain largely 
untested. 

A newly observed phenomenon exhibiting genuine three-dimensional character- 
istics now offers a promising way to probe the coronal dynamics of reconnection: 
dark voids, shown in Figure 3. arc observed by TRACE and Yohkoh SXT descending 
through the haze of 15MK post-eruption loops in the corona (see, e.g., McKenzie 
& Hudson 1999; Gallagher ct al. 2002; Inncs et al. 2003; Asai et al. 2004; Shccley 
et al. 2004). These appear to be truly dynamic phenomena, as they remain coherent 
during their evolution, with dominant, trackable features. In addition, these descend- 
ing voids are clearly three-dimensional phenomena, as their appearance breaks up 
the apparent two-dimensional symmetry of the post-eruption arcade of hot loops in 
Figure 3. A study by Sheeley et al. (2004) even provides evidence that the tracks 
of these voids through the high corona map directly to the tracks of bright loops 
which subsequently appear lower down in the corona. The implication, therefore, is 
that these descending voids are actually evacuated three-dimensional magnetic loops 
which have not yet filled with hot plasma to become visible in emission (McKenzie 
& Hudson 1999). 

These descending coronal voids resemble the teardrop-shaped slow-mode shocks 
from the two-dimensional model of Biernat et al. (1987), making non-steady recon- 
nection a promising explanation of this phenomenon. But the fact that the observed 
voids break the two-dimensional symmetry and then form into fully three-dimensional 
loops means that they must be a three-dimensional manifestation of this model. We 
therefore hypothesize that these descending voids are flux tubes which were gener- 
ated by a short duration, three-dimensional patch of reconnection in the post-CME 
current sheet. Their easily trackable dynamics provide a golden opportunity for mod- 
eling. Our goal is to learn what we can about these coronal voids by extending the 
non-steady reconnection theory into three dimensions. 



- 5 - 



In this manuscript, we focus, as in the corresponding two-dimensional steady- 
state model of Petschek (1964) and the two-dimensional non-steady model of Biernat 
et al. (1987), on the dynamics of reconnected fields rather than the mechanism which 
causes reconnection itself. Sections 2 and 3 explore a simplified scenario whereby 
reconnection and post reconnection dynamics occur in distinct phases. In §2 the 
first phase of brief, high-resistivity reconnection is shown to create a small pocket of 
potential field in the otherwise undisturbed current sheet. This forms two flux tubes 
bent across the current sheet, which then serve as initial conditions for the second, 
post-reconnection, dynamic phase. The dynamical relaxation of these flux tubes in 
the presence of unreconnected flux is developed in §3. This development is done flrst 
in terms of an ideal thin flux tube model and then modified to include additional 
external effects such as "added mass" . Having developed these theories, we then turn 
to three-dimensional magnetohydrodynamical (MHD) simulations to study the details 
of these dynamics in more detail in §4. First, in §4.1, we initialize a set of simulations 
with the analytically reconnected state derived in §2. We study how the retraction 
of the reconnected flux in these simulations compares with the theory of §3. We then 
proceed, in §4.2, to more realistic simulations in which we impose reconnection in a 
one-dimensional current sheet by increasing the resistivity in a sphere on the sheet 
for a short time. In §5, we analyze the added mass effect on the simulated flux tubes, 
estimating the importance of this effect. Finally in §6 we summarize our results and 
briefly discuss their implications for explaining the form and dynamics of post-flare 
coronal voids and arcades. 



2. THE MAGNETIC RECONNECTION EPISODE 

For this study, we assume that the current sheet in which the reconnection occurs 
is one-dimensional, leaving more complex configurations for later studies. For the 
analytical derivation of the reconnected state, we assume the initial magnetic field on 
cither side of the current sheet is uniform, with magnitude Bq, and that the current 
sheet is formed by a tangential discontinuity in the magnetic field at the z = plane. 
The discontinuity has half-angle (, making the explicit form of the field 

B = y-6oCOsC + xSoSgn(2;) sinC, (1) 

where sgn(2;) = ±1 is discontinuous at 2; = 0. This is an equilibrium, provided the 
resistivity is exactly zero. 

As a first step in our study we will consider a hypothetical scenario where the 
reconnection itself and the post-reconnection dynamics occur in distinct phases. That 
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is to say a magnetic reconneciton episode (MRE) occurs in which the resistivity is 
greatly enhanced within a sphere of radius S over a very brief period. To accomplish 
the desired separation we take the period to be so brief that the plasma remains 
stationary both inside and outside the sphere. The resistivity is so greatly enhanced 
that at the completion of the episode the field inside the sphere is completely current 
free: B = Vx. The potential % is harmonic within the sphere, and satisfies Neumann 
conditions to match the radial component of the field of Equation (1). The general 
solution in a sphere centered at Xq = {xo,yo,zo) — can be written in terms of a 
normalized potential Xo(x) 

X(x) = BocosiQy + i?osin(C)5 Xo (|) ■ (2) 

The normalized potential is harmonic within the unit sphere and matches the radial 
component of B/Sq = sgn(2;)x at the outer boundary (i.e., it solves the ( — 7r/2 
case) . 



2.1. Two-Dimensional Case 

To match previous two-dimensional reconnection studies (see, e.g., Petschek 1964; 
Biernat et al. 1987; Nitta et al. 2001; Hirose et al. 2001; Birn et al. 2001; Huba 2005) 
we first consider the case with translational symmetry along y. With this symmetry, 
the reconnection region is a cylinder with axis at x — 0. Writing 

By ^ Bo cos C, (3) 

x5^ + zB, = Bo sin(C)5 Vxo = Bo sin(C)5 y x VAo, (4) 
we can define the complex function 

f(w^^] =Ao-tXo. (5) 



The Cauchy-Riemann equation gives the magnetic field components from the deriva- 
tive of F: 

B, + iB^ = -dy,F. (6) 

To find F, we expand xo as an infinite series of cylindrical harmonics in r = 
(x^ -I- z^y/'^ and (j) — taixr^{z/x), and set its radial derivative at the surface of the 
cylinder (r = 8) equal to the radial component of the magnetic field there: 



dr 



E/j-\m-l 
mAm ( -T ] sin(m0) = sgn(2;)x • f = cos(0) sgn(sin 0). (7) 

r=5 rn=0 
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Using the orthogonality relations for cos(m0), we solve for the constants to give 

m=0 ' 

Using xo = —Im{F), we replace sin(2m0) with w^"* = (e"^r/5)^"^ to find: 

^/ ^ 1 1 - /I + W\ 

F{w) = In , 9 

TT W \1 — W J 

where we have summed over the infinite series. This function is analytic within the 
disk \w\ < 1. 

Contours of the flux function ^0(3:^, z) are shown in Figure 4, giving the projec- 
tions of field lines onto the x — z plane. The field lines form four different flux systems. 
Two of these systems are layers bent back upon themselves around the remaining cur- 
rent sheets. We refer to these two systems as reconnected flux tubes, because they 
consist of field lines which have changed their topology so that they now connect from 
one flux system (at 2; > 0) to a distinct second flux system (at z < 0). On either side 
of these bent layers are regions of flux which remain entirely in a single flux system. 
Some of these fieldines have been modified by the resistivity and, if the guide field 
By is nonzero, have also changed their connections. These could therefore also be 
referred to as reconnected field lines (see, e.g., Hornig 2005), but for this paper, we 
only refer to the field lines which cross between the two flux regions as reconnected 
field lines. 

The reconnected flux per unit length in y in each reconnected flux tube is given 

by the flux crossing the z = plane between a; = and x — 5. This is Ao(w — 
0) - Ao{w = 1) ^ -Re[F(0)] - -Re[F(l)] = 2/7r. Thus the reconnected flux per unit 
length is 

A$ = -BoSsinC- (10) 

TT 

By the same reasoning, the flux entering the reconnection region per unit length is 
Re[F{i)]—Re[F{l)] — 1, or $20 = BqS sin^, and so a fraction 2/7r of the flux entering 
the reconnection region reconnects from one side of the current sheet to the other. 
The remaining flux is distorted by the diffusion, but is not topologically changed, in 
that it docs not cross the current sheet. In this case, the guide field flux does not 
contribute to the total flux entering the reconnection cylinder, as it lies parallel to 
this cyhnder, and so never crosses its surface. 

The current which was formerly within the reconnection region —6 < x < 6, has 
been diffused to the cylinder's surface where it forms a surface current (see Fig. 4). 
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This is not in equilibrium; there is a force at the cyhnder surface, which ultimately 
drives plasma motion. 

The magnetic energy dissipated, per unit length in y, by the resistivity in the 
cylinder is 



= 1 - - / I VxolVdrd^ = 1 - ^ ^ . ^ 0.595 , (11) 

TT J TT^ ^-^ [m'' — 1/4 

m=l ^ ' ' 



Eo,2D 

where Eo^2D = tt^^Bq sin^((^)/87r is the initial reconnection (non guide field) compo- 
nent of the magnetic energy per unit length in y. 



2.2. Three-Dimensional Case 

In three dimensions, the MRE occurs within a sphere of radius 5. Expanding xo 
in spherical harmonics, matching dxo/dr on the surface of the sphere to the radial 
component of the external magnetic field there, and using the orthogonality relations 
for the Legendre Polynomials P„(cos6'), gives 

with r = (a;^ + 1/^ + z'^Y^'^, tan0 = y/x, and cos^^ = z/r. Here P™(cos^^) is the 
associated Legendre function. Using this in the expression for the full field (Eq. 2) 
allows field lines to be traced from a plane y — ya < into and out of the the 
reconnection region, and on to a plane y — yb > S, ets shown in Figure 5. Field lines 
which encounter the reconnection sphere compose two semi- cylinders, on opposite 
sides of ^ = 0. Beginning at y = ya on z > the field lines break into a reconnected 
tube which crosses over to 2; < 0, and a distorted but unreconnected tube which 
remains within z > 0. The separatrix consists of field lines which map to 2; = at 
r = S. 

All field lines in the reconnected flux tube must cross z = within the sphere. 
In the normalized field, the magnetic fiux crossing the half-plane ^ = 0, x > is 



-n/2 

The reconnected flux within the tube which crosses from 2; < to 2; > is therefore 

$i„6e = 5o5'sinC = — sinC, (14) 

TT 
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where $o = BottS'^ is the flux incident on the sphere. The second reconnected flux 
tube, shown in Figure 5, contains the same flux, so 2sin(C)/7r of the flux incident on 
the sphere has reconnected. In the hmit C = 7r/2, this reduces to the two dimensional 
hmit. 

The magnetic energy dissipated by the resistivity is 



where Eq — sin^lQE^S^ /6 is the reconnection component of the initial magnetic 
energy in the sphere. This is very close to the fraction of energy released by dissipation 
in the two-dimensional case. 



Prior to the reconnection episode the magnetic field is in mechanical equilibrium, 
since the magnetic field on both sides of the current sheet is uniform. Following the 
MRE there are two post-reconnection flux tubes with sharp bends; these are clearly 
out of equilibrium. The object of the present work is to model the dynamical evolution 
of these structures, to ascertain how the locahzed reconnection affects the global held. 

Since the reconnection has created flux tubes, one model for their evolution can 
be sought in the dynamical equations for isolated, thin flux tubes proposed by Spruit 
(1981). Such models have traditionally been applied to cases with large plasma P 
where an isolated flux tube is conflned by the pressure of unmagnetized surroundings. 
The underlying approach of considering the net forces acting on segments of a tube 
should, however, be equally applicable to the present situation, where the flux tube 
is surrounded by magnetic field, possibly at very small /3. 

We characterize one flux tube by its axis r, parameterized by its arc-length i. 
The tube encloses a total axial magnetic flux, ^tube, in the form of field with strength 
B{£). The field strength is determined by pressure balance across the tube. In low-/3, 
pressure balance requires that the magnetic field strength within the fiux tube match 
that in the layers between which it is sandwiched: B{£) = Bq. This is the primary 
respect in which the current case differs from the more conventional, high-/5 case. 

A material point, r{i), accelerates due to the net forces acting on it. If the only 
force is the imbalance of magnetic tension, resulting from curvature of the field, the 
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3. 



POST-RECONNECTION EVOLUTION 



3.1. Thin Flux Tube Equations 
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governing equation is 



'dt 



(16) 



where t = dr/d£ is the tangent vector. Writing this in terms of r gives a wave 
equation: 

. (17) 



where v\ = B"^/ (47rp) is the Alfven speed. 



The standard solution to this equation is r = r(£ ± VAt). In general, this is the 
superposition of two waves traveling in opposite directions. As an example, we take 
the shape shown in Figure 6 a, 

r(/,0) =ro-x|£|sinC + y-^cosC, (18) 

as an initial state. The solution to Equation (17) which has this shape and zero 
velocity at t = is 

r(£,t) = ro-x smC + y cosC- (19) 

Taking the derivative of Equation (19) with respect to time gives the velocity 

^VA . ^ f i + VAt i-VAt\ , . 

V = X— sm C -7-^ r + r ■ (20) 

Taking the derivative of Equation (19) with respect to £ gives the tangent vector 

t = -X— ^ — n + ^ ^ smC + ycosC- (21) 



2 \\£ + VAt\ \£-VAt\ 

Together, these two satisfy Equation (16), generating a set of straight segments, shown 
in Figure 6a-c. The central segment moves at velocity v = — xi>^_osinC, while the 
remaining two segments of the tube remain fixed. The corners between segments move 
at a constant speed \d£/dt\ — va,o, where va,q is the Alfven speed in the stationary 
segments. In the case where the initial profile forms the hat shape of Figure 6d, as 
many field lines do in the analytical potential state derived above, the dynamic flux 
tube consist of five straight segments, as shown in Figure 6e-/. 

This analytic solution is a generalization of the two-dimensional solution of Bier- 
nat et al. (1987), but differs in several key respects. The slow-mode shocks corre- 
spond, in this case, to the bends propagating at the Alfven speed. These are not 
simple shocks, owing to the influence of unshocked external field on the tube. It is for 
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that reason that the field strength is unchanged by the bend at \ow-(3. It is a pecuhar- 
ity of the two-dimensional case that makes four different shocks, propagating along 
the four legs of the post-reconnection flux tubes, merge pairwise into two teardrops. 
Lacking this peculiarity, the three-dimensional case has its full compliment of four 
separate shocks. 

Another important feature of this three-dimensional solution is the absence of 
snow-plowing. In the two-dimensional model the retracting flux tube sweeps up mass 
which must be accommodated by an ever-growing bubble (see Fig. 2). In the three- 
dimensional case, however, the bends move apart to create an expanding region in 
which to accommodate the new mass. The moving segment between the bends has 
been shortened by a factor cosC, and so the material on that segment has been 
compressed to a constant density oi p — po/ cos(, while the tube's cross section is 
unchanged (due to the assumed constancy oi B). A version of the two-dimensional 
solution might be recovered by setting ( = 7r/2. However, the divergence in the post- 
shock density is a signal that the general solution does not include the possibility of 
an expanding bubble or cross section. 



3.2. Energy Release 

When the reconnected flux tube retracts, it releases magnetic energy by short- 
ening the length of its field lines. At time t from the rcconncction episode, the shock 
has moved a distance \{t) = v^t along the flux tube in both directions. Or, if the 
shock has reached the boundary by time t, A is simply the original length of the 
flux tube from the boundary to the reconnection site. The segment of the flux tube 
exposed to the shocks by this time has an initial length of 2A, but is shortened by the 
passage of the shocks to 2AcosC (see Fig. 6). This means that the magnetic energy 
per reconnected tube has been decreased by 2$tu6e-6oA(l — cosC)/87r. This energy is 
converted into the kinetic and thermal energy of the moving flux tube segment, in a 
manner analogous to the conversion of magnetic energy by the slow-mode shocks in 
the solution of Biernat et al. (1987). Normalizing energy such that E' = 8ttE/ {Bq^q), 
the energy released by the retraction of both reconnected tubes is 

AE}, = -Asin(C)(l-cosC), (22) 

where we have set ^tube — ^osin(C)/7r (Eq. 14). For a small reconnection region, 
5 <^ X, this is much larger than the energy released in the diffusive reconnection event 
itself: A£;;^^ ~ 0.85sin2(C) (Eq. 15). 
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In addition, the ambient, unreconnected field also releases energy due to the 
reconnection. When the pair of flux tubes retract from the reconnection region, they 
vacate a volume of length 2A(1 — cosC) and area 25^ sin Unreconnected field then 
expands slightly to fill the void, meaning that it must weaken slightly. This external 
field has a pre-expansion energy of 

^ ^ ^o$ext2A(l-cosC) ^23) 

ext StT ■ ^ ' 

Here ^^xt includes all the flux not in the two reconnected tubes, possibly infinite in 
extent: 

^ext ^ Bo{Af^ii-25^ sin C), (24) 

where Af^u is the full area of interest, e.g., the simulation area or a significant fraction 
of the corona. When the reconnected tube retracts, the external field expands to fill 
the vacated volume, decreasing in strength to Bq — b and expanding in area to Afuu 
while keeping the same flux. Setting the initial flux (Eq. 24) equal to the final flux ( 
[Bq — b] A full), we can solve for the change in field strength: 

^ _ Bp A full - Bo{Afuii - 25^ sin C) _ 2BoS^ sin C 

^full -^full 

The energy released by this field is then 

. p 26$e.A(l - cosC) ABo6^ sin (BoiAfuii - 2<52sinC)A(l - cosC) 

^ ^^;:;i ^ ■ ^''^ 

Taking the limit Af^u » 5^, we find that the energy released in the external field is 

A£;;,, = lAsin(C)(l-cosC), (27) 

TT 

which, remarkably, is the same as the energy released in the fiux tube itself as it 
retracts. 



3.3. Inclusion of External Forces 

The foregoing analysis provides an idealized flux tube evolution following recon- 
nection. To obtain the clean analytic form it assumed that the only force acting on 
a tube segment is the magnetic tension of the tube itself. There are in general other 
forces due to the surrounding plasma and un-reconnected field. These forces modify 
the fiux tube dynamics in complicated ways. 
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To estimate these additional effects, we assume tliat tlie solution retains the 
general form of the idealized one: bends propagate at the Alfven speed, creating a 
post-reconnection segment moving with a mean velocity (i.e., center-of-mass velocity) 
of V. The net forces on the post-reconnection segment, or simply "segment", increase 
its momentum. If the segment achieves a steady velocity then 



■net - — ^ V—, (28) 



dP ^dM 
dt dt ' 

where M is the total mass of the segment. As the bends move at speed v\ along the 
two tube legs they sweep up mass at a rate 

^ = 2Po-A^, (29) 

where po is the density of the ambient plasma, and ^ tube/ Bo is the cross sectional 
area of the tube. 

The magnetic tension force along a flux tube, F = $tu6eB/47r, is directed along 
the axis. The contribution in the direction of motion from a single leg is 

-^leg = -^^tubeB±, (30) 

where B± — Bq sin C, is the component of the magnetic field along the direction of 
motion. 

There is also a drag force on the segment, which matches the force the segment 
must exert on the external plasma. It exerts this force in order to deform the un- 
reconnected flux and to accelerate a wake of entrained fluid. We assume that the 
work done by these forces is a multiple X > of the kinetic energy of the segment 
itself. In terms of this factor the net force can be written 

2 dP 

Fnet = -^^tubeB± - (31) 



and Newton's law becomes 



^^tubeB± - 2{X + l)VpoVA^. (32) 
Solving for the velocity of the segment's center of mass gives 

V = = (33) 

4npvA{X + l) X + l ^ ^ 

Thus even as the bends propagate along this direction at fA,±, the center of mass of 
the segment follows behind them at a fraction, 1/(X + 1), of that speed due to the 



-14- 



drag it experiences from its surroundings. This would naturally cause the segment 
to become arched, but otherwise the behavior would resemble that of the idealized 
solution. 

If the principle external force is one of drag, we can estimate the factor X us- 
ing the ratio of kinetic energies of the reconnected flux to the un-reconnected flux. 
Assuming these two components move at similar speed, 

X = (34) 

^rec 

Inserting this expression for X into Equation (33) then gives us the expected velocity 
of the tubes. In §5, we calculate this from our simulations to find out how important 
this added mass effect is. 



4. SIMULATIONS 

We now use three-dimensional MHD simulations to test the predictions derived 
above for the dynamics of isolated magnetic reconnection events. We take these 
simulations in two stages. The first set of simulations reproduce the idealized two- 
stage process of reconnection and relaxation. These use the reconnected state derived 
in §2.2 as the initial condition. The initial current sheet therefore has a sphere 
of potential field on it, as defined by Equation (12). and only a small amount of 
reconnection take place dynamically, due to the background resistivity of Sva/vo ~ 
200. The second set of simulations consider a more realistic case where reconnection 
and relaxation occur together. These start with an unreconnected current sheet and 
impose a sphere of high resistivity, Sva/t} = 2, on the sheet for a short, but finite, 
time. Once the reconnected fiux equals the reconnected fiux in the initial state of the 
potential simulations, this high resistivity region is turned off. 

As these simulations have a finite grid size, we cannot simulate the infinitely 
thin current sheet of Equation (1). Rather, we simulate a finite width Harris current 
sheet, with a guide field component added. The magnetic field of the current sheet, 
shown in Figure 7 a, is then given by 

B^Bq sin(C) tanh (^y j x + -Bq cos(C) tanh (^-j-^ y. (35) 

In this equilibrium, the guide field By drops to zero at the center of the sheet. We 
choose this form to represent the scenario where the two legs of a CME collide to 
form the current sheet, with no fiux initially in between them. This configuration 
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naturally generates a current sheet with a thin, zero field region between the two 
collided flux systems. While we focus primarily on this form of the current sheet, for 
comparison we also simulate two other guide field profiles. The first is where By is 
constant across the current sheet: By = BqCos(. The second is a low-/9 equilibrium, 
where |B| is uniform across the current sheet, requiring that By = {Bq — 5^)^/^. 

In all simulations, we set the current sheet half-width to be the same as the pixel 
size: l/L = 1/64, where the simulation cube is 2L on a side, and the simulations are 
run at 128'^ resolution. The radius of the reconnection sphere is 5/L = 7r/16. The 
pressure is set such that p + B'^/{87i) = constant. Far from the current sheet, the 
pressure is p = po = 20/3 in units where i?o/(87r)^/^ = 1 for the high-/3 simulations, 
and So/(87r)^/^ = 4 for the low-/? simulations. This gives P = Sirpo/B^ = 20/3 and 
5/12, respectively. The density is set to be initially uniform everywhere. 

We use the MHD code CRUNCH3D (see, e.g., Dahlburg & Norton 1995) to 
carry out these simulations. This code solves the compressible, visco-resistive MHD 
equations, in the form presented in Linton & Antiochos (2005), with uniform viscosity, 
but non-uniform resistivity. The code is pseudo-spectral, and therefore has periodic 
boundary conditions. 

4.1. Potential Reconnected Field 
4.1.1. Simulations 

For the first set of simulations, the initial state has a sphere of potential, recon- 
nected field imposed on the current sheet, as shown in Figures 7a and 8a. Figure 7 
is a y = cut through the simulation, showing the magnetic field in the plane as the 
vectors, and the field perpendicular to the plane (the guide field) as the greyscale. 
Figure 7 a shows how the addition of the potential sphere to the current sheet cre- 
ates a set of reconnected fields, which form an x-point in the midplane of the sphere. 
These reconnected fields are well out of equihbrium, as there is no force balancing 
their tension. In Figures 7b and c, these fields therefore quickly pull away from the re- 
connection region, forming a pair of teardrop shapes. These shapes are fully formed 
by Figure 7c, meaning that all of the initially reconnected fiux has retracted from 
the reconnection sphere into the teardrop shapes. From then, the reconnected fields 
simply pull away from the center of the simulation, approximately keeping the same 
teardrop shape throughout. This dynamic resembles the two-dimensional analytic 
dynamics of Biernat et al. (1987) (Fig. 2). 
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There is, however, a fundamental difference between the three-dimensional re- 
connected field shown in Figure 7 and the two-dimensional reconnected field shown 
in Figure 2. While the envelopes of the teardrop shapes are slow-mode shocks in two 
dimensions, this is not the case in three dimensions. To see these shocks, one must 
look at the field lines in a three-dimensional figure, such as Figure 8. Here, a set of 
Lagrangian trace particles are followed dynamically during the simulation, and field 
fines are traced from these particles and plotted at the same times as the slices of 
Figure 7. The trace particles used for this figure initially lie within the reconnected 
sphere along the y — z — axis. These field fines are shown from the — z axis, per- 
pendicular to the view shown in Figure 7. The y — plane of Figure 7 is represented 
in Figure 8a by the white line bisecting the reconnected field lines. The color table in 
this figure is proportional to the field aligned electric current, as shown by the color 
bar. 

Figure 8a shows that, in the initial state, all the field lines from these trace 
particles are reconnected. Most of the field lines take the shape of the upper part of 
a trapezoid, with one bend in the field lines as they enter the potential sphere, and a 
second bend as they exit. The solution for their evolution, given by the thin flux tube 
equations presented in §3.1, should be the superposition of two oppositely propagating 
trapezoids, each of half the magnitude of the initial trapezoid, as in Figures 6d-f. The 
two bends thus propagate in both directions along each field line, giving two pairs of 
shocks on each fiux tube. These are the slow-mode shocks corresponding to that of 
Petschek (1964) and Biernat et al. (1987). One of these pairs is indicated in Figure 
8b-d by the white lines labeled as 'slow shocks.' The field lines do not change their 
shape significantly before the shock hits them, but then they bend to a new direction 
as the shock passes them. These shocks quickly propagate out of the two dimensional 
plane of Figure 7 and therefore do not form the boundary of the teardrop shapes in 
this two-dimensional cut, in contrast to the two-dimensional shocks of Figure 2. 

We performed this simulation, with a potential sphere of reconnection, for angles 
( = [1, 2, 3, 4]7r/8, as summarized in Figures 9a-d, respectively. Here each panel shows 
two snapshots from the evolution of a simulation. The top half of each panel shows the 
initial state of the simulation, and the bottom half of each panel shows the magnetic 
field in the middle of the evolution. The three simulations with a guide field (a-c) 
show the slow-mode shock structure discussed above. Panel d, where the guide field 
is zero, is the closest to the two-dimensional limit, and therefore shows that the field 
lines form a bulging teardrop shape quite similar to the two-dimensional configuration 
predicted by Biernat et al. (1987) (Fig. 2). 
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4- 1-2. Analysis 

According to the thin flux tube theory of §3.1, the reconnected flux tubes should 
retract from the reconnection region at the reconnection Alfven speed va± = va sin (. 
To test this, we measure the speed of the center of the flux tube, i.e., the part lying 
along the y — z — axis. Figure 10a shows a slice of along the y = plane, along 
with a white dashed hne along the y — z — axis. We take the profile of along 
this line at successive timesteps and stack the them vertically to create the distance 
versus time plot of Figure 105. In this plot, features moving in the x direction appear 
as diagonal lines whose slopes give their x velocity. To measure the velocity of the 
center of mass of the flux tube, we flt a gaussian to the x-proflle of positive B^ along 
this line on the x < side of the simulation at each timestep. As an example, the 
centroid of this gaussian is plotted for tvA±/L = 0.658 in Figure 10a as the vertical 
white line. This centroid is then plotted as a function of time in Figure 106 as the 
sohd white hne, while the ±2(7 level of the gaussian are plotted as the two dotted 
white fines. We then find the finear least squares fit to this center-of-mass fine from 
the time when the teardrops are first fully formed {tvA±/L — 0.27) until the fiux 
tube is about to hit the edge of the simulation (tvA±/L = 1.4). The fit, plotted for 
the appropriate time range as the dashed white line in Figure 106, gives a velocity of 
= AOGva^- 

The asterisks in Figure 11a show the measured velocities of the potential recon- 
nection fiux tubes for the four angles simulated. The tube speed is consistently lower 
than the reconnection Alfven speed, Va±, but it approaches this speed as the angle be- 
tween the reconnecting fields increases. In contrast, the slow-mode shock propagates 
at the full Alfven speed, as expected. For ( = n/A, a. distance versus time plot with 
distance along direction of shock propagation, i.e., along the x — y diagonal, gives 
a shock propagation speed of .9blvA- As the thin fiux tube analysis of §2 predicts 
that the tube speed should be vai. if the shock speed is va, it is clear that there is 
a drag force, such as the added mass effect, on the tubes. In addition to the tube 
speed being slower than predicted, a second indication of this drag force is given by 
the fact that the flux tubes of Figure 8 do not form straight segments between the 
shocks, but rather form curved segments. 

Figure 10 shows an additional feature of the three-dimensional reconnection 
which may provide an answer to the source of the drag force: the reconnected tubes 
are sweeping up unreconnected fleld which lies across their path. This field lies on 
or close to the current sheet on z = constant planes, in contrast to the reconnected 
fleld which crosses the ^ = current sheet plane from one set of flux to the other 
set. The two sets of field are therefore topologically entwined, and wrap around each 
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other rather than flowing over each other. The fact that the reconnected field crosses 
the z = plane in one direction forces the unreconnected field it sweeps up to warp 
such that it crosses z = in the other direction. This warped field can be seen in 
Figure 10a. The reconnected field forms the classic teardrop shape, with negative 
on the X < side and positive Bz on the x > side. But in front of these teardrop 
shapes the unreconnected, warped fields are visible as a layer of Bz of the opposite 
sign. 

The second set of simulations wc performed studied the effect of having a stronger 
guide field in the current sheet. In this case, a magnetic equilibrium was simulated 
with (3 = 20/3 as before, but with a uniform guide field, rather than a guide field 
which drops toward zero in the current sheet. If the tubes are slowed down in part 
because they drag the guide field in the current sheet along with them, then the 
flow in this second set of simulations should be slower. Indeed we found that the 
velocities are slower: 34% for C = tt/S, 16% for C = 27r/8, and 6% for ( = 37r/8. Both 
simulations were exactly the same for An/ 8 because at that angle the guide field is 
zero everywhere. Thus, we find that the drag effect becomes stronger as the guide 
field in the current sheet gets stronger. 



4.2. Resistively Reconnected Field 

The potential reconnection simulations of §4.1 allowed us to study the dynamics 
of reconnected fields without worrying about the dynamics of the reconnection itself. 
Having analyzed this simpler situation, we now proceed to study both the reconnec- 
tion and the subsequent dynamics in the same simulation. Here wc start with the 
same unreconnected state given by Equation (35). But instead of inserting the po- 
tential field of Equation (12) into the reconnection sphere, we leave the initial field 
unreconnected, and increase the resistivity in the sphere (r < S) to 

77 = 77o(l + 99e-'-'/'5') (36) 

where rjo is the background resistivity, tiq/Sva = 1/200. The resistivity stays at this 
high level until the amount of dynamically reconnected flux equals that which is 
initially reconnected in the equivalent potential-sphere configuration, and then it is 
decreased back to 770. 
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4-2.1. Two- Dimensional View 

Two-dimensional snapshots of the magnetic field from such a simulation ai ( = 
27r/8 are shown in Figure 12, in the same format as Figure 7. Panel a shows the 
unperturbed current sheet, with contours lines of rj superimposed to show where the 
resistivity is concentrated. As this sheet will be dissipated by the resistive sphere, it 
is interesting to follow its structure on the line x — y — 0, which goes through the 
center of the resistive sphere. The maximum current strength and the width of the 
current at half of its maximum value along this line are plotted as a function of time 
in Figure 13. This shows that the maximum initial value of the current in this sheet 
is 0.7 Jo, where Jq = V x B = Bq sin( is the prescribed initial peak current strength. 
The initial width of the sheet is 2.7 / (where / is one pixel wide). Note that the 
width of the prescribed profile (Eq. 35) is 1.8 / and so the current sheet is wider and 
weaker than prescribed by the initial conditions. This is because the code smooths 
out strongly peaked features with a raised cosine filter as part of the initialization 
routine. 

The effect of the resistive sphere can be seen in Figures 12b and c, where it rapidly 
diffuses the current sheet, causing it to thicken. Figure 13 shows that the sheet width 
jumps to between 12 I and 8 I while the amplitude drops to ~ 0.12 Jq. The diffusion 
causes the magnetic field to reconnect, forming a configuration much like that of the 
potential configuration. This field rapidly pulls away from the reconnection region, 
as before. This can be seen in Figure 12c where the reconnected fiux has started to 
form into two teardrop shapes on either side of the reconnection region. At this time, 
the diffusive sphere is still active, and so new fiux is being continually added to the 
reconnected fiux tubes, and the teardrop shapes extend back to the x-point. 

The fiux reconnected by this diffusion can be estimated by calculating the inte- 
gral of Bz crossing the z = half plane a: > in the positive direction. This is shown 
in Figure 14a as the dash-dotted line. At tvA±/L = .34 as much fiux has reconnected 
(0.8 ^tube) as is reconnected in the initial condition of the equivalent potential re- 
connection simulation. This is smaller than the value, ^tube, given by the analytical 
analysis because the simulated current sheet has finite width I, and so the initial fiux 
$0 incident on the reconnection sphere is only 80% of that expected for an infinitely 
thin current sheet. At this time, we turn the resistive sphere off. Reconnection slows 
significantly at this point, but does not stop, due to the low level of background 
resistivity. 

After the resistive sphere is turned off, the teardrops of reconnected flux continue 
to move away from the reconnection site, and the current sheet reforms behind them. 
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Interestingly, Figure 13 shows that the width and amphtude of the current sheet 
return to approximately their initial values. 

For comparison with this finite duration magnetic reconnection episode, we sim- 
ulated the same interaction but never turned off the resistive sphere. The result, at 
the same time as in Figure 12/ is shown in Figure 125^. In this case, the x-point still 

exists at a; = 0, and the teardrop shape extends all the way back to this x-point. This 
shows that the aspect ratio of the flux tube cross section is strongly influenced by 
the duration of the reconnection episode which creates it. The leading edge of this 
flux tube in panel g has traveled slightly further in the same amount of time than the 
flux tube with a limited amount of reconnection in panel /. This indicates that the 
accumulation of extra reconnected flux in the flux tube increases its speed slightly, 
possibly by helping to overcome the drag force. The current sheet characteristics for 
this steady reconnection simulation are plotted as the dotted lines in Figure 13. This 
shows that the width of the sheet settles at about 8 I and the amplitude stays at 
about 0.12 Jo as the reconnection approaches a steady-state. The reconnection rate 
stays the constant through the rest of the simulation rather than slowing down as it 
does when the resistivity is turned off. 

The flux tube velocities for these high-/9 resistive simulations at the four angles 
C = [1, 2, 3, 4]7r/8 are plotted as the triangles in Figure 11a. These flux tubes are 
slower than the flux tubes retracting from the potential reconnection region at low 
angles, yet they are slightly faster at the higher angles. 

Again, we repeated this set of simulations for a flat guide fleld and found the 
same trend as we did for the potential case: the tubes are slowed by the additional 
guide fleld, particularly as C decreases and the guide fleld gets stronger. Here we 
found that the velocities were 41% slower for ^ = tt/S, 17% slower for C, = 27r/8, and 
12% slower for C = 37r/8. 

4-2.2. Three- Dimensional View 

Figure 15a shows the initial three-dimensional conflguration of this simulation, 
where two sets of straight, unreconnected fleldines cross each other at right angles 
(C = """Z^)- As for Figure 8, these fleld lines are taken from Lagrangian trace particles, 
chosen such that they intersect the sphere of high resistivity at the start of the 
simulation. In this case we show a set which is initially just on either side of the 
current sheet. As the particles lie in the high resistivity region initially, they are 
not frozen onto field lines. They do, however, show snapshots of the state of the 
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reconnecting field, and more exactly follow the field line evolution after they leave 
the resistive sphere, or after the resistive sphere is turned off. 

Initially, the fieldines are perfectly straight and unreconnected. As soon as the 
simulation starts, however, the resistive sphere generates reconnection. The recon- 
necting field lines 'fan' across the current sheet, appearing blue because the current 
sheet has a very high parallel electric current. This fanning agrees with the analy- 
sis and simulations of Pontin et al. (2004) and Pontin et al. (2005). According to 
this theory, in three dimensions reconnecting field lines continuously change their 
connectivity while they are actively reconnecting, rather than undergoing a single 
connectivity- changing event as in two-dimensional reconnection. 

The resulting reconnected fiux tubes progress in a similar fashion to that of the 
potentially reconnected field lines. It is difficult to see the slow-mode shocks in this 
field line figure because the shock fronts are much more diffuse due to the finite time it 
takes for the reconnection to occur as compared with the instantaneous reconnection 
in the potential case. However, the shocks can still be followed in diagonal shces of 
Bz along the x — y plane. Again, we find that these shocks travel at close to the full 
Alfven speed (~ O.Sva for C = 27r/8). 

The four simulations of this type are summarized in Figure 9e-h. The top half 
of each panel shows the reconnected fields early in their evolution, while the bottom 
half shows the tubes in mid-evolution. The tube dynamics are qualitatively similar to 
those of the potential sphere simulations shown in Figure 9a-d. A weak shock front 
can be seen, particularly in the top half of panels e-g, propagating along the flux 
tubes as in the potential reconnection simulations. In addition, the arched shapes of 
the reconnected field lines are qualitatively the same in the two types of simulations. 

4.3. Low-/3 Simulations 

The third set of resistive reconnection simulations we performed are the low-/3 
simulations (/3 = 5/12) where the guide field increases in the current sheet to provide 
pressure balance there. We were unable to simulate the potential simulation at all 
angles at the same level of viscosity and resistivity as for the high-/? simulations, and 
so we do not discuss this experiment here. The problem is that the fields are very 
much out of equilibrium at the start of the potential simulation, and there is less 
gas pressure to balance any sudden disturbances at low- (3 than at high-/?. Therefore 
the density tends to become negative in cavitation regions of strong motion, causing 
the code to crash. We could solve this problem through the artificial method of 
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forcing the density to remain larger than some small value. However, as the resistive 
reconnection simulations do not have the same initial non-equilibrium problem, we 
focus on those simulations instead. 

The two-dimensional view of the low- (3, resistive reconnection simulation at C = 
27r/8 is shown in Figure 16. This shows that decreasing /3 from 20/3 to 5/12 does 
not dramatically affect the reconnection dynamics. The current sheet thickens as 
before when the resistive sphere is turned on, and two teardrop shaped flux tube 
cross sections are formed. The current sheet reforms again once the resistive sphere 
is turned off after Figure 16 c. The velocities of these flux tubes are lower than in 
the high-/5 simulations, as shown by the squares in Figure 11a. In addition, the 
tube velocities are not significantly faster at high ( angles, in contrast to the high- 
P simulations. Finally, for reference. Figure 16g again shows a second simulation 
wherein the resistive sphere is never turned off, and so the fiux tube continues to 
grow in time. This tube again travels slightly farther during the same time than the 
tube of Figure 16/, which contains less reconnected flux. 



5. ADDED MASS EFFECT 

The velocities we derive for all of these simulations pose the question: why are the 
tubes moving at a slower speed than predicted by the thin flux tube model presented 
above? A likely possibility is the added mass effect discussed in §3.3. When a body 
accelerates through an external fluid, as the reconnected flux tubes do here, the 
external fluid is necessarily also accelerated. This takes extra energy, and therefore 
slows down the flux tube. If this is correct, we should see that the unreconnected 
fleld in the simulation gains kinetic energy in proportion to the kinetic energy of the 
reconnected field. 

To study this, we must separate unreconnected fields from reconnected flelds. 
For each time step at which we want to calculate these kinetic energies, we first trace 
a volume filling set of field lines, starting from a uniformly distributed grid of 64^ 
points at the y — plane. To ensure that the field lines are volume filling, we trace 
them in both directions in this plane until they hit the y — ±L boundaries. As By 
is greater than zero at all points in the simulation (for ^ 7^ 7r/2), this routine counts 
each field line only once, as no field line doubles back in the y direction. To fill 
the whole simulation volume with field lines, this requires that we trace field lines 
out of the simulation box in the x direction, making use of the periodic boundary 
conditions. The ( = 7i/2 case is more involved to study, as it has no organizing guide 
field. In this case, we trace all field lines from z > starting at a; = —L (the left hand 
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boundary of the simulation). We trace these field fines until they reach the x — 
plane or the z = plane. Assuming symmetry of By and and anti-symmetry of B^ 
about a; = 0, and symmetry of By and B^ and anti-symmetry of B^ about z = (the 
simulation shows both are valid) we can recover the whole volume of the simulation 
without double counting any field lines. 

Each field line i from the trace is assigned a flux d^i, given by the field strength 
Bi at its starting point (xj, 0, Zi) multiplied by an area 4L'^/N, where is the total 
number of field lines traced. By integrating local quantities along each field line, we 
can then find the field line integrated global quantities, such as volume or kinetic 
energy. The field strength varies along a field line, so we account for the changing 
area along the field line by assuming a constant flux: dAi{x.) = (i$j/|B(x)|. The 
volume of each held fine is then 



As a check of our method, we sum over all N field lines to find the total volume. 
We find that this agrees with the true volume to within .008% for N — 64^ field fines. 
We can then calculate the kinetic energy of each field fine: 



Summing over the all field lines, we find the kinetic energy agrees with the true sum 
to within 3.5% on average. We then must divide the field fines into reconnected 
and unreconncctcd field lines. Here wc take advantage of the simple geometry of 
our simulation: in the unreconnected state, none of the field lines cross the z = 
boundary between the differently directed sets of field lines. Once a field line 
reconnects, however, it connects fields on both sides of this plane, and so it must 
intersect this plane. Thus, to categorize field lines as reconnected, we require that they 
cross the z — plane in the negative direction for x < and in the positive direction 
for X > 0. This rejects field lines which have wrapped around the reconnected tube 
but have not reconnected, as they cross the z — plane in the opposite direction (as 
discussed in §4.1.2). 

The resulting calculations of kinetic energy for the reconnected versus the unre- 
connected fields are shown in Figure 14 as a function of time for the high-/? resistive 
reconnection simulation at ( — 27r/8 (see Figs. 12 and 15). Here one can see that 
the kinetic energy of the reconnected tubes increases in time, and that the kinetic 
energy of the unreconnected field increases in proportion. Thus, the kinetic energy, 
and therefore the velocity of the reconnected tubes must be less than it would in the 
absence of an external medium. 




(37) 




(38) 
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As a check of our field line sorting algorithm, the dotted hne in Figure 14 a shows 
the reconnected flux per tube for this simulation, calculated as half the sum of the 
flux of all reconnected field lines in the field line tracing calculation. This compares 
reasonably well with the dash-dotted line, which shows the reconnected fiux measured 
by summing all the flux which crosses the z = 0, x > plane in the positive direction. 

Wc find that the unreconnected kinetic energy approximately follows the recon- 
nected kinetic energy as in Figure 14b for all simulations studied. To estimate the 
added mass effect, we calculate the ratio of the integrated kinetic energy for the 
unreconnected and reconnected field lines: 

X = ll^fl^, (39) 

where the integrals are taken from the time the flux tubes are fully formed until they 
are about to reach the boundary. We flnd that the ratio X increases as ( increases. 
This is further evidence that the guide fleld is acting as a drag on the flux tubes, as 
the strength of the guide fleld decreases as the angle increases. 

Using the ratio X in Equations 33 and 34, we calculate the expected velocities of 
the reconnected tubes due to the added mass effect. The results are shown in Figure 
lib. Comparing these velocities with the measured velocities of Figure 11a, we find 
that, on average, the potential reconnection tubes move at 82% of this predicted 
speed, the high-/? resistive reconnection tubes move at 72% of this predicted speed, 
but the low-/5 resistive reconnection tubes move at only 49% of this speed. Thus we 
conclude that the added mass effect accounts for most of the flux tube drag at high-/3, 
but there is an additional, signiflcant source of drag at low-/3. 



6. CONCLUSIONS 

We have studied non-steady three-dimensional magnetic reconnection by impos- 
ing reconnection in a localized sphere on a one- dimensional current sheet. We have 
shown that external magnetic fleld intersecting this sphere of reconnection forms 
a pair of bent, reconnected flux tubes. We have solved analytically for the post- 
reconnection evolution of these flux tubes, showing that slow-mode shocks propagate 
as isolated fronts along the legs of the tubes, releasing magnetic energy and acceler- 
ating the fiux tubes away from the reconnection region. These isolated shock fronts 
are the three-dimensional generalization of the two V-shaped shocks of steady two 
dimensional Petschek reconnection (Petschek 1964), and of the two teardrop shaped 
shocks of non-steady two-dimensional reconnection (Biernat et al. 1987). 
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We then studied the evolution of this reconnected magnetic field with three- 
dimensional magnetohydrodynamical simulations. We found that, as predicted, the 
reconnected field quickly retracts from the reconncction sphere, accelerating as the 
slow-mode shock fronts pass by. This accelerated field forms a pair of three-dimensional, 
arched fiux tubes whose cross sections have a distinct teardrop shape. We found that 
the velocities of these fiux tubes are smaller than the reconncction Alfven speed pre- 
dicted by the theory, indicating that some drag force is slowing them down. We 
provide evidence that an added mass effect, wherein the flux tubes sweep up external 
magnetic field and plasma, is largely responsible for this drag force. 

These three-dimensional reconnected fiux tubes are a promising candidate for ex- 
plaining observations of descending coronal voids seen below post-CME fiares. Their 
teardrop shaped cross section is similar to the shape of these voids (McKenzie & 
Hudson 1999), and their three-dimensional arched shape is similar to the shape of 
the coronal arcade loops formed by these voids later in their evolution (Sheeley et al. 
2004). We therefore argue that these coronal voids are fiux tubes formed by localized 
patches of reconncction higher in the corona. If this is the case, then by studying the 
dynamics of these voids, we can directly study the dynamics of reconnected coronal 
fields. To further explore this possibility, we plan to simulate this patchy reconncction 
in a two-dimensional Y-type current sheet. In this case, we expect that the recon- 
nected tubes will decelerate as they descend toward the Y-type null and the arcade 
of field fines below it, much as the coronal voids are observed to do (Sheeley et al. 
2004). We will compare the deceleration of the simulated tubes with the observed 
deceleration of the coronal voids, and therefore more accurately test the validity of 
this model. We will also simulate multiple, simultaneous reconncction events in the 
same current sheet studying how these tubes interact with each other, and whether 
their final equilibrium resembles the tangled, sheared equilibrium formed by post-fiare 
coronal loop arcades. 

We wish to thank Drs. J. Klimchuk, D. McKenzie, N. Sheeley, and H. Warren 
for useful discussions. This work was carried out with support from NASA and ONR, 
and was in part completed while the authors were participating in programs at the 
Kavli Institute for Theoretical Physics and at the Newton Institute. 
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Fig. 1. — Snapshots from a two dimensional MHD simulation of a coronal mass 
ejection, showing the geometry of the fields during eruption. The field bows out 
from the low corona in the left panel, erupting from the solar atmosphere. Once the 
eruption is fully developed, the field starts to pinch off and reconnect on the back side 
of the CME, as shown in the right hand panel. The reconnected fields on the right 
hand side of this newly formed current sheet spring up to the CME, while those on 
the left side spring down into the corona to form arcade loops. From MacNeice et al. 
(2004). 
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Fig. 2. — Theoretical field configuration for a burst of reconnection in 2D. The red 
lines show the slow shock fronts, the green lines show the current sheet, the blue lines 
are reconnected field lines, the black lines are unreconnected field lines, and the pink 
lines are the separatrices between the two. From Longcope (2004), after Biernat et al. 
(1987). 
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Fig. 3. — TRACE Images from the April 21, 2002 X flare. The first three panels show 
unsubtracted images from the flare. The diffuse, bright emission corresponds to high 
temperature plasma from the Fe XXIV 192 A line in the TRACE 195 A bandpass. 
The loop-like structures on the bottom of the arcade are loops that have cooled down 
to about 1 MK and are emitting in Fe XII 195 A. The arrow points to a void that 
is descending through the hot plasma cloud. The final panel is a difference image 
between successive TRACE exposures. From Sheeley et al. (2004). 
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Fig. 4. — Contours of the normalized flux function Aq{x,z), from the real part of 
Equation 9, showing field lines after a two dimensional reconnection episode. The 
reconnection region is enclosed by a thick sold line, inside which the field is current- 
free. Solid hues show the field hues in the two reconnected flux tubes; dashed hues 
show the unreconnected field lines. 
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Fig. 5. — 3D plot of field lines intersecting the sphere of potential reconnected field. 
The field lines shown all cross the ya < plane above the current sheet {z > 0). When 
the field lines pass through the reconnection sphere they deviate from their straight 
trajectory. Some are unreconnected, and so cross to |/b > plane above the current 
sheet, but others have reconnected and so cross the yb > plane below the current 
sheet. The points at which they intersect the two y — constant planes are shown on 
the left and right sides as x — z projections, along with the separatrix curve which 
separates reconnected from unreconnected field lines in those planes (solid lines) , and 
a dashed curve delineating the edge of the flux which intersects the reconnection 
sphere. 
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Fig. 6. — The evolution of two post-reconnection flux tubes according to the thin flux 

tube dynamics of §3.1. A central segment, between the slow mode shocks formed at 
the bends in the tube, moves vertically down due to magnetic tension, while the rest 
of the tube beyond the slow shocks remains fixed. Panels a-c show the motion for 
a fiux tube initially forming the upper part of a triangle, while panels d-f show the 
motion for a flux tube initially forming the upper part of a trapezoid. 
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Fig. 7. — Slices at |/ = of magnetic field from the potential reconnection simulation 
at C = 27r/8> at times tvA±/L = [0.0,0.2,0.4,0.6,0.8,1.0]. The vectors show the 
in-plane magnetic field, while the grey scale shows the guide field {By), where black 
is zero field and white is positive. Panel a shows the initial configuration, including 
the potential sphere at the center. Panels 6-/ show the subsequent evolution, as the 
reconnected fields pull away from this sphere and form into teardrop shapes. 
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Fig. 8. — Three dimensional view of the magnetic field lines from the potential re- 
connection simulation at ( = 27r/8, at the same times as in Figure 7. Field lines are 
traced from Lagrangian particles initially on the y = plane inside the reconnection 
sphere. The color is set by the level of current density parallel to the magnetic field, 
normalized by the maximum initial current density Jq. The slow mode shocks set up 
by the reconnection are visible in panels b-d. The white vertical line in panel a shows 
the location of the y = cuts of Figure 7. 



Fig. 9. — Field lines of the high-/? reconnection simulations at the four angles sim- 
ulated: ( = [1, 2, 3, 4]7r/8, from left to right. The top halves of panels a-d show the 
potential simulations in their initial state, while the bottom halves of these panels 
show the tubes during their evolution at tvA±/L = [0.3,0.6,0.8,0.8], respectively. 
The top halves of panels e-g show the resistive simulations after the same amount 
of flux has reconnected as was initially reconnected in the potential simulations, at 
tvAj_/L = [0.2, 0.3, 0.5], respectively. Because the flux at ^ = Att/S evolves so quickly 
once reconnected, the top half of panel (h) is shown at an earlier time, tvA±/L = 0.3, 
when the reconnected flux equals only half of that reconnected in the potential state. 
The bottom halves of panels e-h show the reconnected tubes later in their evolution, 
at tvA±/L = [0.6,0.9,1.1,0.8], respectively. Note that panels a-c and e-g show the 
field lines from the — z direction while panels d and h, where the guide field is zero, 
show the field lines from the y direction. 
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Fig. 10. — Height versus time stackplot. Panel a shows a y = shce of during 
the potential reconnection simulation at C = 27r/8. A gaussian is fit to the negative 
(black) part of Bz along the white dashed curve for a; < 0, and this is used to find 
the center of mass (shown as the white vertical line) of the teardrop shape. Panel b 
shows cuts along this same dashed line at successive times, stacked on top of each 
other to give a time-distance plot. The time-distance location of the center of mass is 
shown as the solid line in this panel, while the ±2(j level of the gaussian fit is shown 
as the two dotted lines. The linear least squares fit to the center of mass location, 
taken from the time the teardrop shape is first fully formed until it first reaches the 
simulation boundary, is shown as the dashed line. The slope of this line gives the 
velocity of the reconnected flux tube, v — 0A06va±- 
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Fig. 11. — Panel a: measured velocities of the reconnected flux tube velocities. Panel 
b: predictions of this velocity based on the added mass calculation of §3.3. The 
asterisks show the velocities from the high /3 potential reconnection simulations with a 
dip in the guide field at the current sheet. The open triangles show the velocities from 
the high-/3 resistive reconnection simulations with the same dipped guide field. The 
open squares show the velocities from the low-/? resistive reconnection simulations. 
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Fig. 12. — Two dimensional magnetic field slices for the high-/? resistive simulation 
at C = 27r/8, in the format of Figure 7. Panels a through g show slices from this 
simulation at Ival/L = [0.0,0.05,0.34,0.65,0.96,1.27]. This shows the evolution of 

magnetic field which reconnects in a small, 3D sphere in the center of a ID current 
sheet. The contours of resistivity are shown in panel a at 10, 30, 60, and 90 times 
the background level of ^7o/(^'^a) = 1/200. The resistive sphere is turned off after 
panel c, as can be seen by the reformation of the thin current sheet at the center in 
panel d. The 3D view of this simulation is shown in Figure 15. Panel g shows a cut 
at tvA±/L — 1.29 from a different simulation with the same initial conditions, but in 
which the resistive sphere is never turned off. 
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Fig. 13. — Time variation of the current sheet along the x = y = line through the 
center of the high-/3 resistive reconnection simulation of Figures 12 and 15. Panel a 
shows the full width at half of the maximum current value as the solid line. Panel 
b shows peak current strength in the y direction (sohd). The dotted lines in both 
panels show the width and amplitude of the current sheet if the resistivity in the 
sphere is never turned off. 




Fig. 14. — Panel a shows the measured reconnected flux in the high-/? resistive re- 
connection simulation at C = 27r/8, normalized by ^tube from Equation 14. The 
dash-dotted line shows the flux calculated from summing over negative on the 
z — plane at x < 0. The dotted line shows the flux calculated from the held 
line summing routine of §5. Panel b shows the global kinetic energies of reconnected 
(solid) versus unreconnected (dashed) field lines, taken from the fieldlinc summing 
routine. These arc normalized by the magnetic energy of the two unreconnected flux 
tubes initially intersecting the reconnection sphere. 
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Fig. 15. — Three dimensional view of field lines from the high-/? resistive reconnection 
simulation at C = 7r/4, at times tvA±/L = [0.00, 0.30, 0.60, 0.91, 1.22, 1.53]. This shows 
a representative set of field lines which reconnect in the small sphere of high resistivity 
on the initially ID current sheet to form two arched loops. The white line in panel a 
shows the location of the 2D plane of Figure 12. 
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Fig. 16. — Two dimensional view of the magnetic field in the low-/? resistive recon- 
ncction simulation aX ( = tt/A, in the format of Figure 7. Panels a through / show 
this simulation at times tvA±/L = [0.00,0.06,0.34,0.62,0.91,1.20]. In this case, the 
current sheet appears white because the guide field is strongest there (to balance the 
drop in the reconnection component of the field). This contrasts with Figure 7 and 
12, where the guide field drops to zero in the current sheet. The field evolves much 
as it does in the equivalent high-/3 resistive reconnection simulation shown in Figure 
12. Again, the resistive sphere is turned off after panel c, and panel g shows a second 
simulation, in which the resistivity is never turned off, at tvA±/L — 1.20. 



